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Abstract 

Following on from our previous study of the geodesic flow on three dimensional 
ellipsoid with equal middle semi-axes, here we study the remaining cases: Ellip- 
soids with two sets of equal semi-axes with SO (2) x S0{2) symmetry, ellipsoids 
with equal larger or smaller semi-axes with SO (2) symmetry, and ellipsoids with 
three semi-axes coinciding with 5*0 (3) symmetry. All of these cases are Liouville- 
integrable, and reduction of the symmetry leads to singular reduced systems on 
lower-dimensional ellipsoids. The critical values of the energy-momentum maps 
and their singular fibers are completely classified. In the cases with SO (2) sym- 
metry there are corank 1 degenerate critical points; all other critical points are 
non-degenreate. We show that in the case with SO{2) x SO{2) symmetry three 
global action variables exist and the image of the energy surface under the energy- 
momentum map is a convex polyhedron. The case with 50(3) symmetry is non- 
commutatively integrable, and we show that the fibers over regular points of the 
energy-casimir map are bundles over S"^. 



1 Introduction 

The geodesic flow on the ellipsoid is the classical example of a non-trivial separable and 
thus Liouville integrable Hamiltonian system. It is the prime example in Jacobi's famous 
"Vorlesungen iiber Dynamik" [7] and may be considered as his motivation to develop 
Hamilton- Jacobi theory and the solution of the Abel- Jacobi inversion problem. Its modern 
treatment was pioneered by the Zurich school, namely by Moser [TT] and Knorrer [SIE], 
generalising to the n-ellipsoid and providing smooth integrals and the general solution in 
terms of 6'-functions for the generic case of an n-ellipsoid with pair-wise distinct semi-axes. 
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Zung [15] carried out an excellent general approach to the topology of Stackel systems, 
including the geodesic flow on the ellipsoid with distinct semi-axes. In [2] we extended his 
results to the degenerate case. We reviewed the geodesic flow on the three dimensional 
ellipsoid with distinct semi-axes, and analysed what happens in the degenerate case where 
we set the two middle semi-axes equal. We found that the topology of the critical values 
in the image of the energy-momentum map changes and the torus-bundle over the regular 
values is non-trivial and has monodromy. This showed that by making the system simpler 
(i.e. more symmetric) it becomes more complicated (i.e. have a non-trivial torus bundle). 

In this study we complete the classiflcation of the geodesic flow on the remaining three 
dimensional ellipsoids with equal semi-axes, namely the case with two sets of equal semi- 
axes and 5*0(2) x 5*0(2) symmetry, the case with either largest or smallest equal semi-axes 
and 50(2) symmetry, and the case with three equal semi-axes and 50(3) symmetry. In 
section 2 we give a brief review of the geodesic flow on general ellipsoids in order to flx our 
notation. The case with 50(2) x 5*0(2) symmetry is described in section 3. Symmetry 
reduction leads to a (singular) system on the ellipse with an effective potential. The 
symmetry is a torus action which gives rise to a global momentum map. The image of 
the energy surface under this momentum map is a convex polyhedron. Equivalently the 
image of the energy surface under the energy-momentum map is that convex polyhedron, 
except that now the generic flbers are three-tori (by Liouville- Arnold) . The convexity of 
the image is related to results of Atiyah ^ and Guillemin- Sternberg [6] on the convexity of 
the image of the momentum map of torus actions. However, their results apply to compact 
symplectic manifolds only. Generalizations to the non-compact case exist [131 CH] , but our 
example seems to be new. In the end of section 3 we prove the existence of three smooth 
global action variables for this situation. The fact that the third action is smooth (and 
does not have monodromy) is not obvious, since the natrual action deflned via the cycles 
from the separation of variables is only continuous, but not smooth. 

The ellipsoid with the 50(2) symmetry and equal smallest semi-axies is described in 
section 4. The system reduces to a (singular) system on the 2-ellipsoid with additional 
potential and a billiard wall inserted in a plane not containing the umbilic points. Two 
singular values have critical points of corank 1 that are degenerate, while all other singular 
values are non-degenerate. The topology of the flbres is deduced using Poincare sections. 
Section 5 gives similar results for the case with equal largest semi-axes. 

In the final section the two cases with three equal semi-axes are described using singular 
reduction to one degree of freedom. In addition to studying the energy-momentum map 
in this case another interesting object is the energy- casimir map, since the system is non- 
commutatively integrable (or superintegrable). We show that the fibers of the generic 
point in the image of the energy-casimir map is a 2-torus bundle over the 2-sphere. 

2 Review of the the geodesic flow on generic S-eUipsoids 

The geodesic fiow on the generic 3-ellipsoid with distinct semi-axes was described in more 
detail in [2]. Here we will briefiy review some of the terminology to provide a foundation 
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for the examples described here. The standard form of the S-elhpsoid embedded in 
with coordinates x = {xo,Xi,X2,X2,) and semi-axes ^/ai, for < ao < ai < a2 < 0:3, is 

2222 

^1 = ^ + ^ + ^ + ^-1 = 0. 

ao Oil 0L2 oiz 

For the generic non-degenerate elhpsoid the semi-axes are distinct. The Lagrangian of a 
free particle with mass 1 is L = \[x\ ^ x\^ x\^ x\). A Hamiltonian description can be 
obtained by introducing momenta yi = Xi and enforcing the constraint by replacing the 
standard symplectic structure dx A dy by a Dirac bracket. The Dirac bracket preserves 
the ellipsoid Ci = and its tangent space 

XoVo XiVi X2I/2 , x^y^ 
02 — 1 1 1 — u . 

Ctl O2 0:3 

We define also the following notation 

^0 , ^1 , ^2 , ^3 1 ^C*! dC2 

9 2-^ 



Oq af a\ q;| 2 ^ dxi dyi 

We can generalise the constraints Ci and C2, and the factor D, for a n — 1 ellipsoid £ 
embedded in M". Lifting to the cotangent bundle, we have coordinates x = (xq, . . . , Xn-i) 
and conjugate momenta y = {yo, . . . ,yn-i) for T*S embedded in T*R". For this generic 
case the Dirac bracket with Casimirs Ci and C2 is given by 



Xi,Xk}^^ = 0, {Xi,yk}2n = ^ik- , IZ/i, 2/fcl2n, = F|— ' (1) 



where the subscript 2n indicates the dimension of the embedding space T*M". 

Returning to the three dimensional ellipsoid, the Hamiltonian is H = |(yo+2/i +2/2+I/3) 
and the equations of motion are 

x, = {x,,H}s, yi = {yi,H}8, 2 = 0,1,2,3. (2) 

The vector field generated by H is denoted by Xh. 

The system is Liouville integrable with smooth global integrals (in the generic case of 
distinct semi-axes) first found by Uhlenbeck and Moser [TTj 

F.=y^+ E ^"^^y^ - '''y^^\ . = o,...,3. (3) 

On the symplectic leaf of the Dirac bracket given by Ci = C2 = they are related by 
Y^Fi/ai = and they have pair- wise vanishing bracket [TT] . The integrals Fi are related 
to the Hamiltonian by H = l{Fo + Fi + F2 + F3). 

As a preparation for the cases studied here, consider the bifurcation diagram at constant 
energy for the generic case for the four cases in which the semi- axes nearly coincide, as 
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Figure 1: Bifurcation Diagram of almost degenerate Ellipsoids with = (1/4,1/4 + 
£, 1, 2), (1/4, 1/2, 1/2 + £, 2), (1/4, 1/2, 1, 1 + e), (1/4, 1/4 + e, 1, 1 + e), where e = 0.03. 

shown in figure [H The bifurcation diagrams are constructed by separating the variables 
using an ellipsoidal coordinate system and Hamilton- Jacobi theory, as described in [2J. 
From top left to bottom right the cases are close to equal smallest axes, equal middle 
axes, equal largest axes, and equal smallest and largest axes, also denoted by 211, 121, 
112, and 22. For the 22 case we see that the image of the momentum map has only a 
single chamber. For the 112 and 211 cases two chambers are present in the diagrams. We 
must therefore analyse these cases in further detail. 

3 Geodesic flow for the eUipsoid with S0{2) x 50(2) 
symmetry 

The three-ellipsoid £ with two pairs of equal axes is defined by 

2 2 2 2 

rf^ ry^ /y'-' rr>^ 

_0_|__i_|__2._|__3_ = X (4") 
«! «i a2 (^2 ' 

where < ai < a2- The 50(2) x 5*0(2) symmetry group action is given by 

^!{x,y;ei,e2) = {5:,y), (5) 

where 

X = (xo COS 9i — xi sin 9i,xo sin 9i + xi cos 9i,X2 cos ^2 — 2:3 sin 6*2, X2 sin 6^2 + 0:3 cos 6*2), 
y = {yo cos 9i - yi sin 9i, yo sin 9i + yi cos 9i,y2 cos 6*2 - ys sin 6*2, y2 sin 6^2 + ^3 cos 6*2). 
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\1/ is generated by the momentum map A4 = (Ji, J2) : M ^ M^, where Ji = xoVi — xiyo 
and J2 = X2y3 — Xsy2 are the angular momenta. 

The Dirac bracket is given by ^ with = «i and as = a2- 

3.1 Liouville Integrability 

Theorem 3.1. The geodesic flow on the ellipsoid with two pairs of equal axes is Liouville 
integrable. Constants of motion are the energy H = ^{y^ + Z/i + Z/i + vi)' ^'^'^ angular 
momenta Ji and J2. 

Proof. H, J I and J2 commute with respect to the Dirac bracket ([T]). Moreover they 
are independent almost everywhere since they are polynomials and independent at, e.g., 
x = (v/^,0,0,0),y = (0,0,l,0). □ 

We note that we also have smooth globally defined integrals Gi and 6*2- Gi is defined 
by Gi = Fq + Fi, where the Fi are as in ([3]) for the generic three dimensional ellipsoid, 
and letting ao — > ai. Similarly G2 is defined by G2 = F2 + F3 for 02- Hence 

Gi = Vo+yl^ [{xoy2 - X2yof + (xoi/s - x^yo)^ + (xij/a - X2yif + (X2I/3 - X3y2f] , 

«! — ^2 

G2 = y2+yl^ — [{xoy2 - X2yof + {xoys - x^y^Y + {xiy2 - X2yif + (0:21/3 - xsy2f] ■ 

We have the relations 2H = Gi + G2 and 

Gi G2 Ji Jl ^ 

~'~ 22 
ai a2 aj 0:2 

These integrals will be used later in the proof of the non-degeneracy of the singular points 
of the momentum map. 

The group action \E' has the invariants 

'k\ = xI + xI, TTl = y^ + yl, nl = xoyo + xiyi, t^I = Xoyi - Xiyo, (6) 



nl = xl + xl, 7i^ = y^ + yj, = X2y2 + xsys, = X2ys - Xsy2, (7) 

related by 

- {A? - {<? = 0, ^fnl - (nir - (nir = 0. (8) 

^ is not free, but has fixed points at the origin in the xo,xi,yo,yi plane for 62 = and 
in the X2,X3,y2,y3 plane for 9i = 0, or both. These are the only fixed points, and they 
occur for zero angular momentum only. When Ji = nl = ji 7^ and J2 = vr^ = ^2 7^ 
the fixed points are not in ^2), so excluding jij2 = we can reduce by ^ and 

obtain a regular reduced system. 
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Theorem 3.2. A set of reduced coordinates {^1,^2,111^^2) is defined on the reduced phase 
space j2)/'^ , for jij2 7^ 0, by the formulae 

a = a/ttI, i2=\n, r]i = ^=, r]2 = 



T' , . 



The reduced coordinates satisfy the Dirac bracket in M.'^[C,,ri], i.e. {., .}4 as defined in ([T]). 
The mapping R : M.^[x, y] — > M^[^, rj] is Poisson from with {., .js to with {., .}4 and 
the reduced system has reduced Hamiltonian 

+ + + (9) 

Proof. Define a set of coordinate on M'^[^,?7] as sfiown. Tlie new variables (^1, ^2, '^1, ^72) 
are invariant under tlie group action Tlie Poisson property of the map R, i.e. {foR^go 
R}8 = {f,g}4: o R foUows from direct computation of the basic brackets, e.g. {^1,^2)4 = 
{\/xl + xf, \/ X2 + xl}s = 0, etc. Using the identities 



72 72 

+ Vl ^1 = T 

SI S2 



7^2 = 7^ + , ^2 = 7f + ^2' (10) 



the new Hamihonian is found, which is that of the geodesic flow on the elhpse plus an 
effective potential. In the new variables the two Casimirs are 

^ + ^ = 1, ^ + ^ = 0. (11) 

«! a2 «! 0.2 

□ 

The condition jij2 7^ insures that the transformation to rji and the Hamiltonian are 
smooth. The reduced Casimirs show that the reduced phase space is T*S^. Symplectic 
coordinates (0, pd,) may be chosen as 



^1 = v^cos(0), ^2 = A/a2sin(0) 



_ -^/cnsm{(f))p^ _ y/a^ cos{(f))p^ (12) 

ai sin^(0) + a2 cos^(0) ' ai sin^(0) + 0^2 cos^(0) 

These are symplectic coordiantes as {0,p</,} = 1. Expressing the reduced Hamiltonian in 
these coordinates gives 

2 sin^(0) + 0:2 cos2(0) aicos2(0) a2sin^(0)/ 

The contour plot for the reduced Hamiltonian is shown in figure El 

Notice that > 0, so that only one quarter of the circle, O<0<7r/2is the reduced 
configuration space. When allowing the full circle the quotient by the Z2 x Z2 action that 
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Figure 2: Contour plot of reduced Hamiltonian with ai = 1, 0:2 = 2, ji = 1,^2 = 0.4 



flips the signs of .^j, r/j, z = 1, 2 gives back the quarter circle. Here this action is trivial, 
but for the singular reduction it will become more interesting. 

The reduced Hamiltonian has one degree of freedom, and equilibrium points of the 
reduced system in general correspond to two-tori in full phase space. Since the reduced 
Hamiltonian grows monotonously with the momenta, for fixed energy these equilibrium 
points occur at minima of H. Conversely, the image of the energy surface under the 
momentum map has as a boundary the values of momenta for which equilibria are found 
in the reduced system. 

Theorem 3.3. The image of the energy surface £h = {(x, y) G T*S\H = h} under the 
momentum map is the convex region bounded by the polygon defined by the lines 

' ±-p= = v^. (14) 



Proof. The boundary of the image of the energy surface under the momentum map cor- 
responds to the equilibrium points of the reduced flow. Equating {^j, H} to zero gives 

- t:- — + — = 0' (15) 



Aai \ ai a2 

so that r]i = 0, using the Casimir (fTTl) . Using this, the equation {r]i, H} = reduces to 

Dividing this equation by aijfC,f gives = «iji/(«2j2)- Inserting the ansatz = 

±Sy/aiji with a scaling factor s into the Casimir of the ellipse determines s and thus gives 
the result 

= (17) 

j- ji _|_ i2 
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Jl 



Figure 3: Bifurcation diagram for the ellipsoid with oq 
2.h = l. 



OL\ < a2 = «3 for a\ = 1, 02 



The sign in the ansatz is chosen so that ±jj is positive. Evaluating the Hamiltonian 
at these critical points gives the critical values which are the boundaries of the polygon 
shown in figure El □ 

This result is a generalization of a general result by Atiyah [1] and Guillemin- Sternberg 
[6] on the convexity of the image of the momentum map. The classical result applies 
to compact symplectic manifolds only, but extensions to certain non-compact situations 
exist [131 [TUj . The fibres over the momentum map in the compact case follow from the 
general theory. In the non-compact case treated here the fibres are computed explicitly 
using singular reduction. The singular reduction also works when jij2 = 0. However, the 
reduced phase space is not a smooth manifold in this case. Previously we have seen that 
in the regular case jij2 7^ the reduced phase space is the cotangent bundle over an open 
interval in 0, hence it is diffeomorphic to M^. 

The bifurcation diagram figure [3] can be considered as taking the square root twice of 
the diagram in figure [H bottom right. The line Ji = is where the lines Fq = and 
Fi = collapse, and J2 = is where F2 = and F^ = collapse. In the limit coming 
from the generic case the whole of each line would appear to be critical, since they are 
both on the boundary of the image of the energy-momentum map. Recall that the limit 
of -Fo(ao — «i) as ai — » ao (or Fi{ai — ao)) equals Jf, but not Ji. Similarly J| is the limit 
of ^2(^2 — 0^3) ■ Obviously is singular when Ji = 0, but Ji itself is not, and similarly 
for J|. Thus taking the square root twice of figure [H bottom right gives figure El Upon 
this transition the critical points along the lines Fq = Fi = and F2 = -F3 = disappear. 
The boundary of the convex polyhedron the result of taking the square root twice of the 
quadratic curve in figure [H bottom right arising from the double roots. The multiplicity 
of the regular changes from 4 to 1 for every regular point in the image. 
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3.2 Singular Reduction 



Theorem 3.4. The singular reduced space has one or two conical singularities. The only 
singular values of the energy-momentum map for constant energy H = h are the bound- 
aries of the convex polygon forming the image of the energy surface under the momentum 
map. 

Proof. We use invariant theory to perform singular reduction in order to find the reduced 
phase space. The Casimirs for the system expressed in terms of the invariants are 

Ci = ^ + ^- l = 0, C2 = ^ + ^ = 0, (18) 

and expressing the Hamihonian in terms of the invariants gives 

E = \ + vr,^) . (19) 

To find the reduced phase space A^~^(ji, j2)/^ we eliminate the invariants from the 
second relation (IHl). The reduced phase space is then defined by the equations 

a„l /„1n2 -2 „ /i "^^ „2 '^2/„l\2 -2 



together with the inequalities 

7ll > 0, 7tI > 0, TT^ > 0. (21) 

Differentiating the five equations with respect to the six variables and computing the rank 
we find that the surfaces defined by these equations are smooth manifolds unless one or 
both of the ji are zero. When ji = then the non- negativity of irl and tt^ implies that all 
7ri = 0, /c = l,2,3. 

We consider the three cases which may occur. For jij2 7^ we have the strict inequalities 
nl > 0, > O5 ^2 > 0- Thus we can eliminate vrj and the resulting equation is that of a 
two-sheeted hyperboloid. Only one sheet is relevant due to the positivity of the invariants. 
This is the found before using symplectic coordinates (0, p^) . 

For the singular cases it is best to consider a covering of the reduced phase space. In 
[2] this was done for the action generated by a single momentum. Here we repeat the 
construction with both momenta. The result is that we describe the reduced phase space 
as that of the geodesic fiow on the ellipse with coordinates (^,77) quotient by the discrete 
group Z2 X Z2. The group action of the ith factor is given by hipping the signs of ^i,rii. 
When = r/j = is not in the reduced phase space this action merely identifies the four 
quarters of the circle. In the singular case, however, conical singularities appear whenever 
ji = and hence = T^j = is possible. When only one angular momentum is zero this 
gives one conical singularity in the reduced phase space, when both are zero it gives two. 
Thus the reduced phase space for ji = j2 = is the canoe [1] . 
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We consider the energy surface in the reduced phase space. Fixing the energy H = h'm. 
the Hamiltonian f|T9l) . we may use it to ehminate a further invariant in fl20|) because the 
Hamiltonian is hnear in the invariants. Changing h leaves the energy surface invariant up 
to a scahng of momenta, which is always the case for geodesic flows. Fixing one positive 
h gives the equations 

Qi = r^lr^l - {T^lf -jI = 0, Q2 = a, (l - ^\ {2h - vr^^) -^i^Y - j| = . (22) 

V "1/ ' 

The singular points are where VQi x VQ2 = 0. The only real solution to this equation 

is 

-kI = 0, nl = K 23 

Finding the remaining invariants from the Casimirs and energy equation and substituting 
into the relations gives 



2h{7Tl 



TT 



jf = ^^, 3^ = 2ha,[l-^ . (24) 

Eliminating vrj from these equations gives us the equation (fT^ for the convex polygon 
found earlier. In other words the only singular values are on the boundary of the polygon. 
The important thing to note here is that the coordinates axes in the image of the momen- 
tum map are not singular values, unless at the corners of the polygon. In particular also 
the origin ji = ^2 = is not a singular value. The reason is that for h > the conical 
singularities of the reduced phase space are outside the energy surface. 

The non-critical energy surfaces in the reduced phase space are just circles, which can 
best be seen from the Hamiltonian fll3p . In the critical cases they are simply points. □ 

The above theorem confirms that inside the image there are no critical points at all. 
The momentum map can be considered as a slice through the energy-momentum map 
for fixed energy. The complete image of the energy-momentum map is a cone over the 
polygon. Since the geodesic flow does not change with energy (up to scaling) the energy 
can be fixed to some positive value. In the following we are interested in the fibres of the 
energy-momentum map (for that fixed energy), not just the fibres of the momentum map. 



3.3 The Liouville Foliation 

The fibres over a regular value of the energy momentum map is a T^. The fibres of 
singular values are described by the following theorem. 

Theorem 3.5. Points on the edges of the convex polygon are corank one elliptic singular 
points and have have fibres T^. The corners correspond to corank two elliptic- elliptic 
singular points and have fibres S^. 
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Proof. Consider the corner of the polygon given by Jf = 2aih, J2 = 0. Here the invariants 
have the values nl = ai, = 2h, vTg = 0, vr^ = 0, 7r| = 0, 7r| = 0. Take a point in the 
pre-image of the reduction map R from full phase space to reduced phase space 

R : T*R' ^ (vrl, 7cl vr^, vr^, nj, 7r|) (25) 

e.g. X = (v^, 0, 0, 0), 2/ = (0, V2h, 0, 0). For the 50(2) x SO (2) group action the orbit is 

^{x,y;e,,e2) = {x,y) (26) 

where 

5 = (v^oi cos6'i, sin6'i, 0, 0), y = (— V^2/i sin ^^i, -\/2/icos 6'i, 0, 0) (27) 

which is a S*^. Similarly for each of the other corners the fibre is a 5*^. 

The fibres over points on the edge of the polygon are found in a similar manner. Let 
(ji; J2) be a point in the bifurcation diagram on the edge of the polygon, namely a solution 
of ffT^ . but not a corner. Then values of nl, tcI and vTg are given by fpHl) and I^M), and 
the other invariants then follow from the Casimirs and energy equation. In the reduced 
space the invariants for a point (ji, 72) on the boundary are 

7r}=«i(l- Ji), nl = 2h{l~^0=), nl = .[^\j2l = .[^{j^l (28) 

and = 0, 7r| = 0. A point in the pre-image of the reduction map is then found, for 
example, 

x = (y^,0,0,^), i/=(0,^,y^,0). (29) 

Note that as the point (ji, j2) is not a corner then none of the values of Xi, x^, yi or 7/3 
can be zero. Thus the group action of 50(2) x 50(2) on this point gives as an orbit. 

We show the non-degeneracy of the critical points as follows. Consider first of all the 
corners of the convex polygon. We choose the corner for which Ji = 0, J2 = 2a2h to 
show non-degeneracy, but an analogous argument holds for the other corners. For this 
corner we have the critical points Xq = Xi = yo = yi = 0. Here we find Gi = 0, G2 = 2h. 
Evaluating the gradients we have VJi = and VGi = 0, hence the critical points are 
of corank two. The Jacobians of the flow for Ji and Gi are then calculated, i.e. DXj^ 
and DXgi respectively. We evaluate the flows for a single point given by X2 = 0, X3 = 
^/a2, ?/3 = 0, ?/2 = V2h^ and find ^DXj^ + vDXg^- This spans the Cartan sub-algebra; the 
4 eigenvalues are ±i (/x ± 2v^2ha\l (^2 ~ cti))- Hence the critical points are of elliptic- 
elliptic type. 

For the edges of the convex polyhedron, excluding the corners, we proceed as follows. 
Consider first of all the edge for which J\ > 0, J2 > 0. We will give the proof of non- 
degeneracy for this case, but analogous arguments hold for the other edges. Consider the 
new integral K which we define by 
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This is an integral because Ji, J2 and H are integrals. Notice that K = is just the 
equation for the edge of the polyhedron. First we parameterise the edge, for example 
choose 

(ji, J2) = {V2h^iCos\ij), y^2h^2sm\ij)), (30) 
for an angle ip G [0, 7r/2]. 

Then we find a point in the pre-image of the reduction map, e.g. 

= ycii cos(-?/'), Xi = 0, X2 = y/a2 sm{ip) , X3 = 0, 



yo = 0, yi = V2hcos{i)), 1/2 = 0, = V2h sin^ip) . 

Assume further that ip ^ kir /2 ioi k E otherwise the point will correspond to a corner 
of the polyhedron. Calculating the flow of K we get Xk = 0. We are on the boundary 
of the image of the momentum map, the equation of which is = 0. Here we have a 
corank one critical point. The Jacobian of the flow DX^ has two non-zero eigenvalues 

±i2\/2/{a2 + «! + (0:2 — m — ))■ These eigenvalues never vanish and are of 

elliptic type. 

We could equally use general theory [3J to find the singular fibres on the boundary of 
the convex polyhedron, as for a corank r singular point of elliptic type the singular fibre 
is T^-\ □ 



3.4 Actions 

Since the system is Liouville integrable, we may consider action angle variables in the 
neighbourhood of any regular point. Two of the action variables are the angular momenta 
Ji and J2. The third action is the action of the reduced system with one degree of freedom 
obtained from the Hamiltonian (fT3l) along the contours C as seen in Fig. O Hence we 
have the following theorem 

Theorem 3.6. Action Integral. The action variables for the geodesic flow on the three 
dimensional ellipsoid with 5*0(2) x 5*0(2) symmetry are the angular momenta Ji, J2 and 
a third integral I defined by 

I=^i P<^d(l) (31) 

where 



c 



(2h ^-i—- ^:\—] {a, sin2(0) + ^2 cos2(0)) • (32) 

\ q;iCOS^((/)) a2sm (0)y ^ 
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Figure 4: Energy surface in action space with ai = 2,a2 = 2,h = 1. 




Figure 5: Smooth energy surface in action space with ai = 2, 0^2 = 2, = 1. 

The action / is "natural" in the technical sense of [H], and simply means that we 
consider a cycle C such that d9i = d62 = 0, where 9i is the angle coordinate conjugate to 
the angular momentum Jj. Since the action is natural it is not necessarily smooth. 

The energy surface in action space is shown in figure HI Notice that the action / is not 
differentiable with respect to Ji and J2 on the axes ji = and j2 = 0. We will prove 
this by calculating the derivatives of the actions and looking what happens on the axes. 
Further note that the faces of the "pyramid" are not planes. We can smooth the energy 
surface by plotting instead / + | Ji | + | J2 1 . This is shown in figure [5l 

The action integral has simple poles at = 0, ivr with residues ijk, respectively. It 
has branch points at values of for which the momentum vanishes. The integral is 
put into standard form by the substitution s = ^/alCos{(p). The part of the integrand 
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involving the angular momenta is the square root of the following expression 



4 Qij2 - (^231 - 2haia2 2 31^1 



2ha2 



2h 



where G factorises as 



(33) 



(34) 



The polynomial G has si and as roots and hence ±Sj are branch points (in s coordinates) 
of the action integral. In coordinates we have branch points at (pi given by ai cos^(0i) = 



The action / is a complete elliptic integral with Legendre normal form 
/ 



(«? + (a2 - «i) s^) / ^ (fc) - \\ ^T^rfen (a^ A;) 



+ 



{si - ai) {al + {a2 - ai) sf) 



ai {al + {a2 - ai) sj 



2\1/ 



17^^^(/5^A:) 



(35) 



where S and 11 are Legendre's complete elliptic integrals of the second and third kind 
respectively and the parameters are 



p _ ("2 - ai){sl - si) ^2 



al + {a2 - ai)sl 



a 



al 



aia2 



-k\ (36) 



s\{a2-ai} {ai - sl){a2 - ai] 

The derivatives of the action / with respect to the other action variables Ji and J2 are 



dl 



dl 



■h 



2-Kai 



\ 



ai sin^(0) + a2 cos2(0) 



2h 



■-'1 



72 

■^2 



Qi cos^ (0) 02 sin ((/)) 



COS^(0) 



J2 



27ia2 



\ 



ai sin^(0) + a2 cos'^{(f)) dcj) 
2}i 5 3 sin^fi 



sin^(0) 



A, 



P2. 



(37) 



(38) 



dl 



has a simple pole at = | whereas ^ has a simple pole at = 0. A picture of the 
complex plane with the branch points and poles is shown in figure [6l For Ji = we have 

ai and < < ai. In other 



Sl 



and < S2 < «!, and for J2 



we have s^ 



words, when Ji = the branch point at s = Si coincides with the simple pole at s = 
(i.e at = ^) and when J2 = the branch point at s = S2 coincides with the simple pole 
at s = a^^ (i.e. at = 0). The integration paths in the complex plane are deformed so 
that the integral may be split up into two separate integrals for each case. The integral 
Ci around the branch points is expanded into a loop Bi around the pole in question and 
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a) 



+ 



+ 



+1 





B2 



57 



Figure 6: Complex plane C(s) and choice of branch cuts (a), integration paths for pe for 
— > (b) and s| — > cti (c). (d) and (e) show decomposition of Ci for case (b) and C2 
for case (c). 
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then the contribution from the pole Ri is subtracted. For the case where si coincides with 
we have 



+ f (39) 



1 /2 

and for the case where S2 coincides with we have 



(40) 



Evaluating the residues of the integrand at the simple pole in question, we have 



Res/3i = — ^ and Resp2 = -77^^- (41) 
Zm I J 1 1 Zm I J 2 1 

We then have 

dl , Ji , , dl ^ J2 
lim — — = lim - — 7 and lim — — = lim — r (42) 

Ji^o dJi Ji^o\Ji\ J2^o dJ2 J2^o IJ2I 

by the method of residues. Hence 

dl , dl .91 , 9/ , , 

= 1' ^™ ^ = ^ = ^™ ^ = 1 • (43) 

Ji-»o+ aJi Ji-»o- aJi J2-»o+ 0J2 J2-»o- aJ2 

Denote the natural actions in each of the four quadrants of the plane in (Ji, J2) as /++, 
/_+, / and /+_. They are related by 



I — \-{~'Jly J2) — 'S'i/++(Ji, J2), 
/__(-Jl,-J2) = 52/-+(-Jl,J2), 
/+_(Jl,-J2) = 5i/__(-Ji,-J2) 
-^++(•^1) •^2) = >S'2/++(Jl, — J2), 



(44) 





where 

^1=1 1 I , ^2 = I -1 I . (45) 

\ooi/ \ooiy 

We then define unimodular matrices Mj such that the appropriate natural actions and 
products of actions and Mj join smoothly at the Ji and J2 axes. We have 



(46) 



I++ 


= Mi7_+ 




= Mi7++, 


Jl 


= 0, J2 


>0, 


h- 


= M2/__ 


= M25l/+_ 


= M2/+_, 


Jl 


= 0,J2 


<0, 


I~+ 


= Ms/— 


= M3-52/-+ 


= M3/-+, 


J2 


= 0, Jl 


<0, 






= M4,52/++ 


= M4/++, 


J2 


= 0,Ji 


>0, 
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on the appropriate axes. 

We see that Mi and M2 have an eigenvalue 1 with eigenvector (0, J2, /)* and M3 and 
M4 have an eigenvalue 1 with eigenvector (Ji,0,/)*. The eigenvector equations and the 
fact that Mj e SL{3,Z) show that Mj must have the form 

/100\ / 1 0\ 

Mi=\ m 1 \ ,ie{l,2}, 1 Lie {3,4}. (47) 

V A 1 y V A 1 / 

In fact, we have Mi = M2 and M3 = M4. As the derivatives of the actions must join 
smoothly at the axes we have 

which implies that 

Ml = M2 = 1 I , M3 = M4 = I 1 I . (49) 





Taking a loop around around the origin and calculating the overall monodromy matrix 
M, gives 

M = (M4,S2)-'(M25i)-'(M3>52)(Mi5i) = id, (50) 

so monodromy is not present in the system. Globally smooth action variables can be 
defined by /, Mi^iJ, M^S2MiSiI and {M2Siy^ M^S2MiSiI in each of the four quadrants 
of the plane, where I is the natural action in the (+, +) plane. Thus we have proved 

Theorem 3.7. The geodesic flow on the three dimensional ellipsoid with two pairs of 
equal axes has three global actions. 

Two of these actions are trivial, they are the generators of the symmetry. The third 
action is given by a non-trivial complete elliptic integral, and when using the separating 
coordinate system this action is not globally smooth. Instead it is only continuous across 
the hues = 0. Nevertheless, even though this natural action is not differentiable, 
different linear combinations exist which give actions that arc smooth across the lines 
Jk = 0. Moreover, completing a cycle around the origin there is no overall monodromy. 
Therefore a globally smooth and single valued third action exists in this case. It would 
be interesting to check whether it is possible to find a single formula in terms of standard 
complete elliptic integrals which is smooth. 

The absence of monodromy was to be expected since there are no critical points inside 
the image of the momentum maps. Nevertheless, from the point of view of the explicit 
calculation of the third momentum this is a non-trivial observation. 
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4 Ellipsoid with 112 symmetry 



Consider the geodesic flow on a three elhpsoid with the two largest semi-axes equal, namely 
< ao < ai < «2 = ^3, corresponding to the symmetry group action S0{1) x S0{1) x 
S0{2). The action of S0{1) has no effect as it merely fixes the variables on which it acts, 
but we include it in a direct product to label the case where we have two equal semi-axes 
so as to distinguish whether it is the two largest, the middle two or the two smallest 
semi- axes which are equal. We will call the 112 ellipsoid the three dimensional oblate 
ellipsoid, in analogy to the two dimensional ellipsoid of revolution with a S0{1) x SO (2) 
symmetry, which is oblate or disc shaped. The Casimirs Ci and C2, the resulting Dirac 
bracket ([I]) and the Hamiltonian are the same as in the generic case, with 02 and as set 
equal. In this situation the integrals F2 and F3 from the generic case are not defined any 
more, but the singular terms cancel in the sum G = F2 + F^. The other integrals Fq and 
Fi remain the same. The system is invariant under rotations in the (x2, ^3) plane and its 
cotangent lift, the 50(2) group action being 

^x,y-9) = {x,y) (51) 

where 

X = {xq, xi, X2 COS 9 — X3 sin 9, X2 sin 9 + X3 cos 9) ^^^^ 
y = iVo, 1/1, Z/2 cos 9-y3 sin 9, y2 sin 9 + y^ cos 9) . 

The group action $ is the flow generated by the angular momentum J = X2yz — X3y2, 
which is a global action variable. 



4.1 Liouville Integrability 

Theorem 4.1. Liouville Integrability. The Geodesic flow on the ellipsoid with equal 
largest semi-axes is Liouville integrable. Gonstants of motion are the energy, H = 
I (1/0 + ?/i + 1/2 + ^/i); angular momentum, J = X2I/3 — X3y2, and the third integral 
G = F2 -\- F^ 

2,2, (2^o?/2 - X2yof {xoys - x^yof {xiy2 - X2yif {xiy^ - x^yif 
<^ = 2/2 + 1/3 H \ \ \ (53) 

0-2 ~ O.Q (y.2 — tto 0-2 — Oil Cy.2 — ttl 

Proof. The relationship between the integrals and the Hamiltonian is 2H = Fq + Fi + G, 
and in the limit as — * 0:2 the constant of motion (02 — a^)F2 becomes J^. Hence G and 
H commute because the Fi commute in the generic case, and commutes with H and 
G, and hence so does J. Here the relationship between the constants of motion on the 
symplectic leaves of the Dirac bracket is given by 

^ + - + --4 = 0, (54) 

ttO «1 Ci2 Ct2 
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which is the hmit of the generic relation YlPi/^^i = using ^2/0^2 + -^3/03 = G/a^ — 
^2(^2 - a3)/(a2tt3)- 

H, J, G and the Casimirs Ci, C2 are functionally independent almost everywhere on the 
level set Ci = C2 = 0: They are polynomial and independent e.g. at x = (^00,0,0,0), 
l/=(0,0,v^,0). □ 

The group action $ has the invariants 

TTi = X2 + X3, 'K2 = yl + yl, vTg = X2y2 + x^y^, 7r4 = a;2?/3 - a;3?/2 , (55) 

related by vri7r2 — TTg — 7r| = 0. The remaining variables Xq, Xi, t/q, |/i are trivial invariants of 
$. The only fixed point of $ is X2 = X3 = 1/2 = I/3 = 0. As before, when J = 7r4 = j 7^ 
this fixed point is not in J~^{j) and the reduction by the S0{2) symmetry leads to a 
smooth reduced system on J^^{j)/ S0{2): 

Lemma 4.2. A set of reduced coordinates (6, ■Ci? 6; ^70) ^71) "'72) is defined on the reduced 
phase space Pj = J^^{j)/ S0{2) by the formulae 

^0 = xo, 6 = xi, G = v^, ^70 = I/O, Vi = 2/1, ^72 = 



The reduced coordinates satisfy the Dirac bracket in M.^[C,,t]], i.e. {., .}q as defined in ([T]). 
The mapping R : R^[a:, ?/] — > M^[6?7] is Poisson from with {., .}§ to with {., .jg and 
the reduced system has reduced Hamiltonian 



1 f 



and additional integral 

^2 ~ C^O ^2 ~ ^2 V Ol2 — tto ^^2 " CH 

Proof. Define a set of coordinates on M^[6 r^] as shown. The Poisson property of the map 
i?, i.e. {f o R,go = {/, g^^ o R follows from direct computation of the basic brackets, 
e-g- {6,6}6 = {2^1, a/x^ + xl}s = 0, etc. The reduced Hamiltonian H and additional 
integral G are found by writing them in terms of the invariants and then expressing them 
in terms of the reduced variables. □ 

The reduced system is the "geodesic flow" on the 2-dimensional ellipsoid with semi-axes 
^/aQ, a/oi, ^/a2 and an additional effective potential j^/2^|, and as ^ > 0, the reduced 
system for |j| > has only the open half of the ellipsoid as configuration space with the 
plane ^ = dynamically not accessible. 

We define a singular coordinate system on M^[xj,?/j] by 



2^2 = 6 cos ^, X3 = 6sin6' (56) 
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where 9 is the angle of rotation corresponding to the <S'0(2) symmetry group action $. 
We then choose confocal elhpsoidal coordinates and parameterise the 2-elhpsoid to give 
us a set of generahsed coordinates (Ai,A2,^^) on the three cUipsoid with the larger two 
semi-axes equal. We denote the conjugate momenta arc denoted by {pi,P2,P9)- 

Lemma 4.3. The Hamiltonian for the geodesic flow on the ellipsoid with equal largest 
semi- axes in local symplectic coordinates reads 

^ _ 2(q;o - Ai)(q;i - Ai)(a2 - Ai) 2 2(ao - A2)(ai - A2)(Qi2 - A2) 2 
Ai(A2 — Ai) A2(Ai — A2) 

_^ («o - a2)(ai - 0-2) 2 
2q;2(Ai - q;2)(A2 - a2)^^' 

The constants of motion are pe and Gi 

~ _ 2(q;o - Aj)(ai - Ai)(Q;2 - Ai) 2 (ao - a2)(«i - "2) 2 
— ^ Vi — f^^i ^ — 7T ^ Pe 

where i — 1, 2. The integrals G and G^ are related by 

^^^^^^ {a2-ao){a2-ai)G _ ^^^^ ^ al - apa, 
a2 

Proof. The Hamiltonian is found by applying the cotangent lift to the new coordinates. 
The variables are separated by multiplication by A2 — Ai and rearranging to find the Gi. 

By separation of variables the momenta Pk conjugate to A^ can be expressed as 

where Q is a cubic polynomial. □ 
The relation of the separation constant G to the constant of motion G is 
Fo ^ F, ^ G ^ J' _Q{z) 



z — ao z — ai z — a2 {z — a2)^ ^i^) 



4.2 Singular Reduction 

Lemma 4.4. The singular reduced phase space of the geodesic flow on the 3-ellipsoid with 
equal largest axes is the phase space of the geodesic flow on the 2-ellipsoid reduced by the 
Z2 action S{^o,^i,^2,Vo,Vi^V2) = (■Co, ^1, -^2, Vo,Vii -^2)- Thus it is the geodesic flow on 
the 2-ellipsoid with a hard billiard wall inserted in the ^2 = 0-plane. 
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Proof. Singular reduction for j = leads to a reduced system on a non-smooth manifold. 
We use invariant theory to carry out the singular reduction. The reduced phase space is 
a subset of ]R''[a;o, yo, Xi,yi, ni, tt2, tt^] and is defined by the two Casimirs 

_o ^ + ^ = 1^ + .JUl + _^ = 0, (58) 

ao cti a2 tto tti a2 

the relation between the invariants ttitT2 — 7r| = and the inequalities tti > 0, > 0. 
Since the Casimirs are linear in the invariants we can eliminate them from the relation to 
give 



ai 1 vTa - + = J , (59) 

which defines a four dimensional object. This is the reduced phase space Pj = J^^{j) / S0{2) 
which is a smooth 4 dimensional manifold for j 7^ and not a smooth object for j = 0. 

The geodesic flow on the 2-ellipsoid reduced by the Z2 action S is the billiard. In [2j we 
proved that the singular reduced phase space of the geodesic flow on the 3-ellipsoid with 
equal middle semi-axes and vanishing angular momentum j = is the phase space of the 
geodesic flow on the 2-ellipsoid reduced by a Z2 action. For the 112 situation the result is 
analogous, although here the Z2 action fixes the plane orthogonal to the largest axis, i.e. 
.^2 = 0, whereas in the 121 case the Z2 action fixed the plane orthogonal to the middle 
axis, namely .^1 = 0. The reader is referred to [2] for a full description of the proof. □ 



4.3 The Liouville foliation 

We investigate the topology of the invariant level sets obtained by fixing the constants 
of motion. The energy momentum map is SA4 = {H, J, G) : M ^ R^. Since H for a 
geodesic flow is homogeneous in the momenta we can fix the energy to, say, h. 

Theorem 4.5. The image of the energy momentum map EM. for constant energy H = h 
is the region in bounded by the quadratic curves (see figure^ 

9 = h -] , (60) 

0:2 - «! a2(«2-tti) 

2^2 , "0 .2 (aA\ 

h 7 rJ , (61) 



a2 — ao a2(a2 — ao 



8a2h C(2 ~ c^oc^i 



(a2 - aoj(a2 - aij a2(a2 - aoj(a2 - aij 

Singular values of the energy momentum map are the boundary curves (elliptic), their 
transverse intersections (elliptic- elliptic), the boundary between the two chambers of the 
bifurcation diagram (hyperbolic) , and the non-transverse intersections between the bound- 
ary curves (degenerate). 
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Proof. As in the generic case critical points can occur on sub-ellipsoids. On xq = yo = 
the integral Fq = and V-Fq = 0, similarly for xi = yi = and Fi. In both cases the 
corresponding sub-ellipsoids are ellipsoids of revolution. The image of the critical points 
with = I/O = is found using the relation flMl) to ehminate Fi in 2H = Fq + Fi + G, 
which gives 

2H = G-aA — -^ 

and hence the upper boundary curve of critical values flUUl) . A similar computation for 
critical points with xi = yi = gives the lower curve flUTj) . which divides the bifurcation 
diagram into two chambers. 

Since for Xi = yi = the integrals F^ = and also its gradient vanishes, G = F2 + F^ 
and its gradient clearly vanishes when X2 = X3 = y2 = ys = 0. Considering the Casimirs 
the solutions set of 0:2 = = ?/2 = 2/3 = is a geodesic flow on the ellipse in the xq-xi 
plane. Fixing the energy two critical circles are obtained. On these critical points also 
J = so that the origin in the image (j, g) = (0, 0) is a critical value. This is analogous 
to the case with equal middle semi-axes. However in this case the critical value is not 
isolated but lies on the intersection of the two curves defined by flU2|) . The curves flU^ 
arise in the following way. 

We use the ellipsoidal coordinates from lemma 14.31 to establish which other points are 
non-singular. These coordinates are non-singular outside any sub-ellipsoid = 0. To 
find critical points in the region of phase space with ^ coordinates such that all 7^ 
we compute the rank of the matrix D{Gi, G2,pe), see lemma. 14.31 Since the variables are 
separated this implies Pi = and 2ha2{Xi — a-if = (^2 — ao)(tt2 — C(i)Pg > so here we 
have critical points not contained in the coordinate singularities = 0. Solving for Aj 
and substituting this value and Pi = into the expression for Gi in lemma 14.31 gives the 
following curves in (j, g) coordinates 



~9 = -ha2 ± |^-| / 2M^^-^o)(^2-^i) (63) 

^2 

and using the relationship between g and g in lemma HiB] to convert to {j,g) coordinates 
we get the curves (l62l) . For the ellipsoid with equal middle axes the equation arising here 
had no solution. The curve intersects tangentially with the curve corresponding to critical 
points Xi = yi = at the points 




I 1 2ha2(.a2 ~ C(i) 2h ^ \ 1 

{j,g)=\±\l^^—^,j^—-y,{a,-2aoa2 + aoa,) | (64) 

Now let us investigate the types of critical points. For the corank one critical points 
on the upper boundary curve flUUl) these points are non-degenerate because the Jacobian 
DXfq of the flow generated by Fq restricted to the critical points xq = yo = is given by 
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the following non-zero sub-block with i = (all other sub-blocks having zero entries at 
the critical points) 



Note that KQ{y,y) never vanishes because ~ cto in the denominator of each term is 
always positive. The matrix is traceless so the square of the eigenvalues is given by the 
negative determinant. Evaluating the matrix on the point xq = xi = X3 = 0, yo = 1/2 = 
shows that the eigenvalues of this matrix never vanish and are of elliptic type, and hence 
these critical points are non-degenerate. 

For the corank one critical points on the lower boundary curve f l61l) we have the Jacobian 
of the flow generated by Fi restricted to the critical points Xi = j/i = as given by flU^ . 
However here Ki{y,y) can vanish because the terms in the sum can have different signs. 
If we consider the point xq = = on all critical sets we have X2 = ±y^a2 and ?/2 = 0. 
The eigenvalues of the matrix block vanish at y^ = 2h{ai — ao)/('^2 ~ c^o) so here we 
have degenerate points. These corresponds to the two points of tangency where the 
lower boundary curve (lUTI) and third curves intersect, defined by ( IMI) . For values of 

< \J^^^^^^^^ we have y^ > 2h{ai — ao)/(tt2 — cto) and here the eigenvalues are real 
and hence hyperbolic. Note that this is the case even for the point on the axis j = 0. For 

> we have y^ < 2h{ai — q;o)/(«2 — «o) and the eigenvalues are purely 

imaginary and hence elliptic. 

The two corank two points given by the intersection of the two curves ( 160|) and (l6Til 
are also non-degenerate, because the non-zero 2x2 blocks of the Jacobians are distinct, 
so that fiDXpg + vDXp^ spans the Cartan subalgebra; the 4 eigenvalues (for any point 
on the critical circles given by Xq = Xi = y^ = yi = 0) are ±2in^2aQh / {a2 — oiq) and 
±2ii/v^2ai7i/(a2 — 0.1). This orbit is a relative equilibrium, i.e. a circle in the X1-X2 
plane. The eigenvalues of DXp^ and DXp^ are elliptic, so at their intersection an orbit of 
elliptic-elliptic type is found. 

Now consider the type of singularities on the third curves flU^ . These are corank one 
critical points. We consider the Jacobian matrix of the flow generated by the one degree 
of freedom system corresponding the the constant of motion G. The eigenvalues A are 
given by = and both partial derivatives are strictly positive. Hence critical 

points on these lines are elliptic. For the origin on these lines we have J = and G = 0. 
This is a critical point of corank two. Moreover VJ = as well, and the Jacobian of 
Xj has eigenvalues ±i, since its flow $ is a rotation. We have that fiDXc + vDXj 
spans the Cartan subalgebra; the 2 eigenvalues (for any point on the critical circles) are 
±i/i 1/80:2/?./ (0:2 — a;o)(a;2 — oti) ± iv. Hence the origin is of elliptic-elliptic type. 

This establishes the existence, non- degeneracy or degeneracy, and type of all critical 
points. The bifurcation diagram is shown in figure [71 

It remains to show that there are no other critical points. We check the pre-images 
in full phase space of the sub-ellipsoids = where the ellipsoidal coordinates are not 




where 




(65) 
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Figure 7: Bifurcation Diagram for the ellipsoid with ao < «i < 02 = as for ao = 1, ai = 
2, q;2 = 3, /i = 1. 

defined. Here the argument is that given for the case of the ellipsoid with equal middle 
semi- axes in [2]. 

□ 

The bifurcation diagram figure [7] can be considered as the square root of the diagram 
in figure [1] bottom left, similar to what happened in the case with two equal middle 
semi- axes [2]. 

We now find the fibres of the energy momentum map at the singular values in the 
bifurcation diagram. 

Theorem 4.6. The singular fibres over the upper boundary curve (!60|) of the image of 
the energy momentum map at constant energy, with the exception of its intersections, 
are two-dimensional tori T^. At each intersection point of the upper and lower boundary 
curves ( |60l) and (I6T1) the singular fibre is . The singular fibre at the origin is two circles 
. For the boundary curves ( |62l) except at its intersections the singular fibre is two sets 
of two-dimensional tori T^. For the lower boundary curve fl6Tl) the singular fibre is a two- 
dimensional torus on the boundary of the bifurcation diagram, B x T"^ on the region 
which divides the bifurcation diagram into two chambers, and for the two degenerate 
points where the curve is tangential to ( !62l) . where B is the Bolsinov-Fomenko atom. The 
multiplicity of the regular fibres T^ is one for the chamber enclosed between the upper 
boundary curve ( 160|) and the lower boundary curve f l6Tl) and two for the chamber enclosed 
between the lower boundary curve (1611) and the curves (lU^ . 

Proof. Consider first of all the upper boundary curve (1601) with Fq = of the image of 
the energy momentum map. On this curve all singularities are of elliptic type, and hence 
the singular fibre is T^~^ where r is the corank of the singularity; hence we have on 
the curve and at its intersection points with the curve (|6Ti) . 

The curve consists of all orbits in the geodesic flow on the ellipsoid of revolution defined 
by Xq = yo = 0, which is an oblate two dimensional ellipsoid, obtained by deleting the 
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smallest axis in the three dimensional ellipsoid. Reduction maps each of this system 
to a relative periodic orbit. The isolated periodic orbit in the 23-plane of the geodesic 
flow on the ellipsoid of revolution corresponds to the extremal points with J = ±y/2a^. 
Reduction maps this relative equilibrium to the fixed point ^ = (0, 0, ^/a2) on the larger- 
axis of the reduced ellipsoid. We also note that the multiplicity of the fibres is one 
on this boundary. If we consider the fibres at the intersection points we find that at 
one intersection we have the periodic orbit in the 23 plane, in full phase space, in one 
direction J = +v^2a^, and at the other the other periodic orbit in the opposite direction 
J = —\p2ooJi^ so the multiplicity is one at the intersections. Moving from the intersection 
to the boundary curve Fq = the multiplicity of must therefore also be one. 

Consider now the other transverse intersection of boundary curves at G = 0, J = 0. 
This is an elliptic-elliptic point of corank two and the fibre is IS^ . The fibres are the 
periodic orbits in the 01 plane on the ellipsoid in full phase space corresponding to vri = 

= = 0. There are two periodic orbits; one in either direction. 

Now consider the curves defined in fl62|) on which the origin (J, G) = (0, 0) lies. For 
points on the curves away from the tangential intersection with flUTj) the critical points are 
of elliptic type and of corank one. Hence the singular fibres are T^. Since the multiplicity 
at the origin is two, the multiplicity on the boundary curves must also be two. 

Finally consider the lower boundary curve ( |6T]) corresponding to Fi = 0. The curve 
consists of all orbits in the geodesic flow on the ellipsoid of revolution defined by xi = yi = 
0. This is an oblate two dimensional ellipsoid, as occurred for the curve corresponding to 
Fq = 0, but in the current instance the ellipsoid of revolution is obtained by deleting the 
second smallest axis in the three dimensional ellipsoid. This is why hyperbolic behaviour 
is possible on this line, and also why the line can have degenerate points. For the part of 
this curve which is a boundary of the bifurcation diagram, excluding the two degenerate 
points, the critical points are elliptic and corank one, and by the same argument as for 
the curve (l60l) the singular fibres are of multiplicity one. 

Before looking at the remaining part of (I6T1) let us consider the multiplicity of the 
regular fibres. The regular fibres are by the Liouville- Arnold theorem. For the chamber 
enclosed between the upper boundary curve flUUl) and the lower boundary curve flUTj) the 
multiplicity is one because the multiplicity is one on the boundaries. For the chamber 
enclosed between the lower boundary curve ( 16T1) and the curves ( l62l) the multiplicity is 
two as the multiplicity is two on the boundaries. 

As we cross the part of the lower boundary curve (|6T1) which divides the bifurcation 
diagram into two chambers, the multiplicity of the regular fibres changes from one to 
two. From the Fomenko theory of atoms we therefore require an atom corresponding 
to a bifurcation from a multiplicity of one to two. A likely candidate for the singular 
fibre in the reduced phase space is the direct product of and the Fomenko atom B. 
We shall show that this is indeed what we get. The critical points are of corank one 
but are of hyperbolic type. In the reduced ellipsoid the critical sets are circles in the 02 
plane. However, the singular fibres also have a separatrix because the critical points are 
hyperbolic. 
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Figure 8: The intersection of the preimage of points on (I6T1) for j 7^ with the Poincare 
section xq = in reduced phase space. The separatrix is of type B. 




Figure 9: The intersection of the preimage of points on (I6T1) for j = with the Poincare 
section xq = in reduced phase space. The separatrix is of type C2. 
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By lemma 14.41 the reduced system for j = is the geodesic flow on the 2-ellipsoid 
quotient by the Z2 action S. Ignoring the quotient the reduced singular fibre consists of 
the unstable periodic orbits corresponding to ^1 = ?7i = on the reduced ellipsoid and 
their separatrix. This can be seen from the Poincare section ^0 = 0. Since ^0 = ''?o = is 
an invariant subflow the boundary of the section with rjo > is an invariant set and it is 
the only place where the flow is not transverse to the section. In configuration space the 
section condition is the ellipse in the 12-plane, and it can be parametrised by an angle 
by (^1,^2) = (yoi COS0, ^/a2sm(j)). The momentum conjugate to then gives the 
momenta as {rji,rj2) = {— ^/al sin (j), y/a2 cos (j))p^/d where d = aisin^(0) +a2cos^(0). 
The reduced Hamiltonian can be solved for rjo on the section and thus the integral G can 

be written as a function of on the section. Setting q = — 7?'^, where 

we assume h = 1, and rearranging the expression to find we have the section of the 
singular fibre is given by 

2 ^ f 2/1^2 f 

^'^ V"2-ao (02 - "i) sin^(0) 

Figure [H] shows the section for j 7^ for > which is the Fomenko atom B. There 
is also a mirror image of B if we plot the part of the graph for which < . However, 
for j = we get a different picture. This is shown in figure [9l Note that here we have 
the Fomenko atom C2. This is consistent because if we start with the graph for which 
j 7^ 0, then as we decrease the value of j the two copies of the atom B approach each 
other, finally merging into one curve at j = 0. If we increase j then as we approach the 
degenerate points the graphs of the B atoms in the section reduce in size until all that is 
left of them at the degenerate points is two points in the (0, p^) plane on the p,p = axis. 

Now the quotient with respect to 5* has to be performed. In the new coordinates the 
Z2 action S is (0,p<^) — > {—(p.—p^). This action has no fixed points on the curve. For 
j' = 0, under the action the C2 is converted into B. To see this, choose a fundamental 
region of the plane, > say, then glue to —p^ at = and similarly at = tt. 
For j 7^ the two B merge into one so we again have B. 

Since the reduced flow is transverse to the section on the singular fibre the complete 
reduced singular fibre is B x S^. The singular fibre in full phase space is found by letting 
$ act on this set. There are no fixed points under $ so every point will be multiplied by 
S^. Hence the singular fibre is 5 x T^. 

For the degenerate points recall that the B atoms reduced to two points which under 
the residual action map to a single point. The singular fibre in reduced phase space is 
therefore S^. In full phase space we have a singular fibre of T^. This is logical because as 
we travel along the boundary curve, where the fibre is T^, we reach the degenerate point 
and the fibre does not change. This is expected behaviour for a degenerate point 
When we enter the interior of the bifurcation diagram and the critical points become 
hyperbolic, the fibre changes to B x T"^. So we always have a in the fibre along the 
boundary curve. 

□ 



\ {a2- ao){a2- ai)d 2/,n 

^ cos . (66) 

/ a2{d-ao) 
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As in the case of the 121 elhpsoid we can define action variables for the 112 situation 
in terms of hyper-elhptic integrals. However, we do not have an isolated singularity in 
this case nor do we have monodromy. The fibre bundle over the regular points in each 
chamber of the bifurcation diagram will be trivial as in the generic case with distinct 
semi- axes. 



5 Ellipsoid with 211 symmetry 

Consider the geodesic fiow on a three ellipsoid with the two smallest semi-axes equal, 
< ao = Q^i < «2 < tts, corresponding to the symmetry group action S0{2) x 5*0(1) x 
5*0(1). We call this ellipsoid the three dimensional prolate ellipsoid, in analogy to the 
two dimensional ellipsoid of revolution with a 5*0(2) x 5*0(1) symmetry. The results are 
analogous to the 112 case, so we will only briefiy summarise them. The integrals Fq and 
Fi from the generic case are not defined any more, but the singular terms cancel in the 
sum G = Fq + Fi. The other integrals F2 and F3 remain the same. The system is invariant 
under rotations in the (xo,Xi) plane and its cotangent lift, the 5*0(2) group action being 

^x,y;e) = {x,y) (67) 

where 

X = (xo COS 9 — Xi sin 9, Xq sin 9 + xi cos 9, X2, X3) ^^^^ 
y = (yo cos 9 - yi sin 9, yo sin 9 + yi cos 9, y2, ys) ■ 

As before the group action $ is the fiow generated by the angular momentum J = 
xoyi — xiyo, which is once again a global action variable. 

Theorem 5.1. Liouville Integrability The Geodesic flow on the ellipsoid with equal small- 
est semi-axes is Liouville integrable. Gonstants of motion are the energy H = ^ (^/q + yf + y^ + vi), 
the angular momentum J = x^yi — xiy^, and the third integral G = Fq + Fi 

n 2,2, (a;oZ/2 - a;2l/o)^ , {.^\y2 - X2yif {xoys - x^yof , {xm - x^yif 

G = yo + yi^ \ \ \ (69) 

Oq — OL2 tto — 0.2 OLq — OL-i tto — OLz 

Proof. The proof is analogous to that of (14. ip . □ 
The group action $ has the invariants 

7l^=xl + xl, 'K2 = yl + yl, TTg = Xo?/o + = Zo?/l - Xi?/o , (70) 

related by 7ri7r2 — 7r| — 7r| = 0. The remaining variables are trivial invariants. The fixed 
point of $ is Xo = Xi = I/O = Z/i = 0. As before, when J = tia = 3 1^ ^ this fixed point 
is not in J~^{j) and the reduction by the 50(2) symmetry leads to a smooth reduced 
system on J-^{]) / S0{2): 



28 



Figure 10: Bifurcation Diagram for the ellipsoid with ao = ai < 02 < as for ao = 1, 0:2 = 
2,a3 = 3,h = 1. 



Lemma 5.2. A set of reduced coordinates {C,o,^i,C,2,Vo,ViyV2) is defined on the reduced 
phase space Pj = J^^{j)/ S0{2) by the formulae 

^0 = VtH, ^1 = X2, 6 = X3, VO = Vl = 1/2, V2 = 2/3- 



The reduced coordinates satisfy the Dirac bracket in M.^[C,,ri], i.e. {., .}q as defined in ([T]). 
The mapping R : M^[x, ?/] MP[^,ri] is Poisson from with {., .jg to with {., .}q and 
the reduced system has reduced Hamiltonian 

1 f 

and additional integral 



{^0V2 - ^2V0? ^ - ^SVO)' + A + ^ ^3^ 

ao — a2 ao — a^ \ ao — a2 a^ — a^ 



Proof. This is analogous to the 112 situation. □ 

We get lemmas analogous to lemma 14.31 for the Hamiltonian in local coordinates and 
lemma H3] for the singular reduced space. For the Liouville foliation we have: 

Theorem 5.3. The image of the energy momentum map EM. for constant energy H = h 
is the region in bounded by the quadratic curves (see figure [W\) 

2ao , , "2 .2 

9 = h + — -J , (71) 

ao — a2 ao(a2 — aoj 

" h+ , , J^ (72) 



ao - as ao(a3 - ao) 
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2 



8aoh a2«3 — Oq 



("2 - tto)(a3 - "o) ao("2 - ao)("3 - "o) 



J- (73) 



Singular values of the energy momentum map are the boundary curves (elliptic), their 
transverse intersections (elliptic- elliptic), the boundary between the two chambers of the 
bifurcation diagram (hyperbolic), and the non-transverse intersections between the bound- 
ary curves (degenerate). 

Proof. This is analogous to the 112 case. □ 

Note that the curvature of the boundary curves in figure [TU] is different to the curvature 
of the boundaries in figure [7] for the 112 elhpsoid because the boundary equations have 
different coefficients. For the 211 elhpsoid, G will always change sign due to the form 
of the quadratic equations for the boundaries, which are always negative for J = and 
always positive for the maximum J value. The tangency points can be positive or negative 
depending on the values of cto, ^2, «3- 

The fibres of the energy momentum map at the singular values in the bifurcation 
diagram are then classified by 

Theorem 5.4. The singular fibres over the boundary curve (ITTj) of the image of the 
energy momentum map at constant energy, with the exception of its intersections, are 
two-dimensional tori T^. At each intersection point of the boundary curves (171 p and (1721) 
the singular fibre is S^. The singular fibre at the origin is two circles . For the boundary 
curves (1731) except at its intersections the singular fibre is two sets of two-dimensional 
tori T^. For the curve f l72|) the singular fibre is a two-dimensional torus for each 
point on the boundary of the bifurcation diagram, B x T^ on the region which divides the 
bifurcation diagram into two chambers, and for the two degenerate points where the 
curve is tangential to ( I73l) . where B is the Fomenko atom. The multiplicity of the regular 
fibres T^ is one for the chamber enclosed between ( I7T|) and ( 1721) and two for the chamber 
enclosed between ( 172|) and ( 1731) . 

Proof. This is analogous to the 112 case. □ 



6 Ellipsoids with 13, 31 symmetry 

Consider the geodesic fiow on the three dimensional ellipsoid corresponding to the group 
actions 50(1) x 50(3) and 50(3) x 50(1). Taking 50(1) x 50(3) first we have 

™2 /v»2 ^2 ^2 

^ + ^ + ^ + ^ = 1 (74) 

tto «1 ttl ttl 

where ao < ai. The Casimirs and the Dirac bracket are the same as before with cti — 
a2 = as- 
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For systems with non- commutative integrability, a more useful concept than the energy 
momentum map, when describing the dynamics of the system, is the energy Casimir 
map SC. This is defined here as £C : T*R^ ^ R x R, (x, y) ^ {H, J) where H = 
I iVo + 2/1 + 2/2 + Vs) is energy and J the total angular momentum arising from the 
symmetry group action of SO (3). 

Theorem 6.1. Non-commutative Integrability. The equations of the geodesic flow on 
the ellipsoid with a symmetry corresponding to 5*0(1) x 5*0(3) are non- commutative in- 
tegrable; namely two independent involutive integrals are the energy H and total angular 
momentum J, and two independent but non-involutive integrals given by angular momenta 
Li2,Li3. The fibre over a regular point of the energy Casimir map {H, J) = {h,j) is a 
fibre bundle with base space 5*^ and fibre T^. 

Proof. The 5*0(3) group action gives rise to three angular momenta L12, L13, L23 (defined 
by L12 = Xiy2 — X2yi etc.) belonging to the Lie algebra so(3)*. These are all functionally 
independent. The total angular momentum J is defined by 

•^^ = -^12 + -^13 + -^23- i"^^) 

Taking J as an integral leaves two functionally independent Lij. By a direct calculation 
{H,J} = 0, {H,L,,} = 0, {J,Lij} = 0, {^12,^^13} = 2:21/3 -a;3l/2. (76) 

Hence we have two independent involutive integrals and two independent but non-involutive 
integrals. We apply Nekhoroshev's theorem [12] on non-commutative integrability to show 
that motion takes place on invariant two tori T^. However, this is when we fix the con- 
stants of motion. Fixing J, L12 and L13 is equivalent to fixing L12, -L13 and L23, assuming 
that they are in the correct range (i.e. obviously L23 is defined by L23 = J"^ — L\2 — L\j^ and 
this must be positive). Now the individual angular momenta Lij belong to so(3)* which 
is isomorphic to R^ with the usual isomorphism of multiplication being mapped to the 
vector product. But this allows the Lij to vary over so is not quite the situation here. 
For the energy Casimir map where we fix (iJ, J) = {h,j), the total angular momentum 
J is fixed so the Lij are free to vary but are constrained on a sphere 5^. So each regular 
fibre of the energy Casimir map is a fibre bundle with fibre and base space 5^. □ 

An analogous situation to this result was described by Fasso [3] for broadly integrable 
Hamiltonian systems, of a meadow of actions in which there are flowers whose petals are 
tori parameterised by angles conjugate to the actions and whose centres are coadjoint 
orbits. 

For the group action 50(3) we have the invariants 

711 = xl + xl + xl, (77) 

7r2 = yl + yl + (78) 
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TTs = xiyi + X2y2 + xsVs- (79) 

The variables Xq, Uq are trivial invariants of the group action. We use the invariants to 
define variables in a reduced system as follows: 

Lemma 6.2. A set of reduced coordinates {C,o, CiyVo^Vi) defined on the reduced phase 
space Pj = J^^{j)/ S0{3) by the formulae 

c c I — ^3 



The reduced coordinates satisfy the Dirac bracket in M^[,^,?7], i.e. {., .}4 as defined in ([T]). 
The mapping R : M^[x, y] — > r/] is Poisson from with {., .}§ to with {., .}4 anc? 
t/ie reduced system has reduced Hamiltonian 



representing a one degree of freedom system on an ellipse. 

Proof. The proof is analogous to lemma I4.2[ □ 
We can now carry out singular reduction by 

Lemma 6.3. The reduced phase space Pj = J^^{j)/ S0{3) is an open subset ofT*S^, 
diffeomorphic to M^, in the case that j ^ 0. For j = the reduced space is two dimensional 
with a conical singularity. 

Proof. The Casimirs and the Hamiltonian are linear in the invariants. The reduced phase 
space is a subset of M^[a;o, yo,TTi,n2, rr^] and defined by the Casimirs, the relation 7ri7r2 — 
^3 = ^"^^ the inequalities tti > 0, > 0. Eliminating the invariants tti and tt^ using 
the Casimirs gives 

ai 1 7r2 = ] , (80j 

which defines the reduced phase space Pj = J^^{j)/ S0{3). As before this is smooth 
except for j = 0. For j 7^ we have > and the reduced phase space is one sheet of 
a two sheeted hyperboloid, which is diffeomorphic to M^. We can also see this by looking 
at the two Casimirs in the reduced variables (Co; Ci; ''?0; which give an open subset of 
T*S^. We only have an open subset and not the whole cotangent bundle as d > 0. For 
j = we have 712 > and here we get a conical singularity in the reduced phase space. □ 

For the singular fibres: 

Theorem 6.4. Singular Fibres For J = we have a singular fibre S"^ x and for 
= 2aih we have a singular fibre S0{3). 
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Proof. For the geodesic flow on the elhpsoid corresponding to a symmetry group 5*0(1) x 
5*0(3) we have < < 2aih. Eliminating the invariants in the relation using the the 
Casimirs and the Hamiltonian gives 



a,{l-^) {2h-y',)-^^i = J' (81) 



Differentiating this relation with respect to the variables xq and yo and setting the deriva- 
tives equal to zero gives a critical point at Xq = 0, = 0. Here the invariants are tti = ai, 
H2 = 2/i and tts = 0. This gives us the singular value = 2aih. Inserting this value into 
(EH) gives 



(fti - Qq) 2\ _ 

J- H Xn — 



and hence the only solution is Xq = 0,yo = 0. Thus the singular fibre is defined by the 
equations 

xl + xl + xl = «! (83) 



yl + yl + yl = 2h (84) 

xm + X2I/2 + a^3Z/3 = 0. (85) 
But this is just the unit circle bundle over 5*^ which is known to be 5*0(3). 

Note that the relation we differentiated to find the singular points was equal to J^. This 
means that another singular point arises due to the problem with differentiability at J = 0. 
Setting (ISTI) equal to zero and rearranging gives 

Xr\ ~\~ ~ 

The curve represented by this equation is homoemorphic to a circle 5*^ in reduced space. 
For any point on this circle we choose any point in the inverse image of the reduction map 
and look at the orbit of 50(3). As the point is never a fixed point of the group action, 
or we would require {xi,X2, X3) = (0, 0, 0) and (yi, 1/2, 2/3) = (0, 0, 0) which is not possible, 
the orbit is always 5*^ and hence the singular fibre is 5^ x 5^. □ 

For the 50(3) x 50(1) group action, the situation is similar to the previous case. We 
have 
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Theorem 6.5. 5*0(3) x 5*0(1) Ellipsoid. For the geodesic flow on the ellipsoid corre- 
sponding to a 50(3) X 5*0(1) symmetry, the behaviour of the flow is analogous to that of 
the 5*0(1) X 50(3) ellipsoid. For total angular momentum J we have < < 2aoh. 
For J = we once again have a singular fibre 5^ x 5*^, for = 2aoh we have a singular 
fibre 5*0(3), and for intermediate values of the points are regular and the regular fibres 
are a bundle over 5*^. 

The ellipsoid with a 5*0(4) symmetry is just the geodesic flow on 5*^, which is well 
known. It has been described in detail by Cushman 
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